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Abstract A bar rotating in a pressure-supported halo generally loses angular mo- 
mentum and slows down due to dynamical friction. Valenzuela & Klypin 
report a counter-example of a bar that rotates in a dense halo with little 
friction for several Gyr, and argue that their result invalidates the claim 
by Debattista & Sellwood that fast bars in real galaxies require a low 
halo density. We show that it is possible for friction to cease for a while 
should the pattern speed of the bar fluctuate upward. The reduced 
friction is due to an anomalous gradient in the phase-space density of 
particles at the principal resonance created by the earlier evolution. The 
result obtained by Valenzuela & Klypin is probably an artifact of their 
adaptive mesh refinement method, but anyway could not persist in a 
real galaxy. The conclusion by Debattista & Sellwood still stands. 
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1. Introduction 

It is now well established that a bar rotating in a halo loses angular 
momentum through dynamical friction. This topic has received a lot of 
attention recently for two important reasons: (1) it offers a constraint 
on the density of the DM halo (Debattista & Sellwood 1998, 2000), and 
(2) it may flatten the density cusp (Weinberg & Katz 2002). 

Both these claims have been challenged. Realistic bars in cuspy ha- 
los produce a mild density decrease at most (HoUey-Bockelmann et al. 
2003) or even a slight increase (Sellwood 2003), but we leave this issue 
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Figure 1. The time evolution of the bar pattern speed in a number of resimulations 
of model Ai of VK03. The evolution reported by VK03 is reproduced as the dot- 

dashod line; all other lines are from simulations from the same initial particle load, 
but run with our code using many different sets of numerical parameters. 



aside here and concentrate instead on the density constraint. Holley- 
Bockelmann & Weinberg (2005) announce a prehminary report of sim- 
ulations with weak friction in halos having uniform density cores, but 
we focus here on the older counter-example claimed by Valenzuela &: 
Klypin (2003; hereafter VK03) of a bar that experiences little friction in 
a cusped dense halo. 

VK03 kindly made available the initial positions and velocities of all 
the particles of their model Ai , in which the bar did not slow for 2-3 Gyrs 
after it had formed and settled. We have used our code (Sellwood 2003) 
to rerun this simulation many times, and the pattern speed evolution 
in many of these runs is shown in Figure 1. It is striking that in most 
cases, the bar slowed earlier than VK03 found, but in one anomalous 
case, the bar stayed fast for about 10 Gyr! The anomalous result is not 
a consequence of some inadequate numerical parameter, since many of 
the other cases are from models with parameters that bracket those of 
the anomalous case - i.e. longer and shorter time steps, coarser and finer 
grids, etc. 

Note that apart from the crucial delay in the onset of friction in the 
case by VK03 and the one anomalous case we find, the evolution is 
generally very similar. In particular, whenever the bar slows, it slows 
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at a similar rate. The following sections account for the discrepancies 
between the results shown in Fig. 1. 



In a classic paper, Tremaine &; Weinberg (1984) laid out the math- 
ematical apparatus for friction in a spherical system. Following the 
precepts of Lynden-Bell &; Kalnajs (1972), they derived a formula for 
the torque experienced by a rotating perturbation potential, <l>p. They 
work in action-angle variables (see Binney &: Tremaine 1987, §3.5). In 
a spherical potential, there are two non-zero actions: the total angular 
momentum per unit mass L = and the radial action J,., each associ- 
ated with two separate frequencies, Q, and /«, which are generalizations 
to orbits of arbitrary eccentricity of the usual frequencies of Lindblad 
epicycles familiar from disk dynamics. In the limit that a constant am- 
plitude perturbation rotates steadily at ilp, they showed that the net 
LBK torque is 



TLBK oc ^ im— + k—\\^rnnk?'5{n9, + kK-rnVlp), (1) 



where / is the usual distribution function and ^mnk is a Fourier coeffi- 
cient of the perturbing potential. The Dirac delta function implies that 
the net torque is the sum of the separate contributions from resonances, 
where nO. + kn = mfip. Because the bar pattern speed decreases, as a 
result of the frictional torque, this expression needs to be generalized to 
a time-dependent forcing (see Weinberg 2004), but the revised expres- 
sion for the torque still contains the same derivatives of the distribution 
function. 

Lynden-Bell (1979) ofi^crcd a clear insight into how an orbit is affected 
when close to a resonance. The unperturbed orbit, which is a rosette in 
an inertial reference frame, closes in any frame that rotates at the rate 



for any pair k, m. [See e.g. Kalnajs (1977) for illustrations of several of 
the most important shapes.] When the pattern speed of the bar is close 
to J7' for some pair k, m, the orbit can be regarded as a closed figure 
that precesses at the slow rate 
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Q' = Q + kK/m, 



(2) 



fig = i'^' - ^p) < ^p- 



(3) 



Under these circumstances, the "fast action" is adiabatically invariant, 
while the "slow action" can suffer a large change. Things are particularly 



Figure 2. Left: The time evolution of the bar pattern speed in the restricted 
simulation discussed in §3. This simulation employs lOM particles. Right: The mean 
density of particles as a function of Lres at t = 200 in the same simulation. 

straightforward at corotation, where the fast action is the radial action, 
while the slow action that can suffer a large change is simply the angular 
momentum. 

3. Restricted Simulations 

Fully self-consistent simulations are complicated by evolution of the 
total potential, changes to the bar mass profile, etc. It is therefore easier 
first to try to understand "restricted" simulations in which a rigid bar 
rotates in halo of non-interacting test particles (Lin Sz Tremaine 1983; 
Sellwood 2004). The particles move in a rigid halo potential that is 
perturbed by that of the rotating bar, and the bar is accelerated in 
response to the vector-sum of the non-axisymmetric forces felt by the 
particles. Figure 2 shows an example of the pattern speed evolution of 
a homogeneous ellipsoid with principal axes 1 : 0.5 : 0.1 in a Hcrnquist 
(1990) halo. The bar mass is 1% of the halo mass, M/j, has a semi-major 
axis equal to the halo break radius r/j, and rotates initially so that it 
just fills its corotation circle. We use units such that G = Mh = r/j = 1. 

At intervals during the simulation, we compute Jl' = + kn/m for 
every particle, and calculate F, the average density of particles as a 
function . It is somewhat easier to understand a plot of F as a function 
of the angular momentum, Lres, of a circular orbit that has the given 

The right-hand panel of Figure 2 shows the form of F near to corota- 
tion (m = 2, = 0) at i = 200, which is typical. Near to the resonance 
(marked by the vertical line), particles cross corotation in both directions 
on horse-shoe orbits (Binney &: Tremaine 1987, §7.5). The generally neg- 
ative slope of F implies an excess of lower L particles that gain angular 
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Figure 3. The solid curve shows evolution of the bar pattern speed in a restricted 
experiment in which the bar was reaccelerated (by external interference) between 
times 180 and 200, but was otherwise allowed to evolve freely. (This experiment used 
IM particles.) The dotted curve shows what happens without interference. 



momentum and move out across the resonance, and this imbalance is 
responsible for friction. If the pattern speed were to stay constant, the 
imbalance would tend to flatten the slope of and the distribution 
of particles about the resonance would approach kinetic equilibrium in 
which there would be more nearly equal numbers of particles crossing 
in both directions. But as fip declines, the resonance keeps moving to 
larger Lrcs, and equilibrium is never established; instead, the density of 
particles about the dominant resonance(s) responsible for friction takes 
on the characteristic humped form shown in Fig. 2. 

Friction arises principally at corotation over most of the evolution. 
The outer Lindblad resonance contributes in the early stages, but domi- 
nates only if the bar is unreasonably fast. The inner Lindblad resonance 
becomes important only when the bar is already quite slow; e.g. it is 
responsible for the more rapid braking around t = 1000 in Fig. 2. 

4. Anomalous Situation 

Now suppose that ftp rises for some reason, after having declined for 
some time, as illustrated in Figure 3. The shoulder in F created by 
the previous friction survives, but the resonance at the now higher Clp 
lies on the other side of the shoulder, as shown in Figure 4. Thus the 
local gradient in F has changed sign, leading to an adverse gradient 
for friction, and the torque from the dominant resonance disappears. 
Under these circumstances, a balance between gainers and losers is soon 
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Figure 4- The mean density of particles as a function of Lres at several different 
times in the simulation shown in Fig. 3. 



established, and the bar can rotate in a dense halo with little friction, 
which wc describe "metastable state". 

In fact, ^Ip declines slowly because of weak friction at other resonances, 
and normal friction resumes when the slope of F at the main resonance 
changes, as shown in the last frame of Fig. 4. 

5. Self-consistent Simulations 

If we now re-examine Fig. 1, we see that the period of weak friction is 
preceded by a small rise in the bar pattern speed in both the simulation 
of VK03 (dot-dashed line) and in the anomalous case we found. It is 
likely therefore that friction stopped for a while in both cases because 
the local density gradient across the principal resonance became flat, as 
just described. 

Analysis of our simulation that displayed this behavior suggests that 
ilp rose because of an interaction between the bar and a spiral in the 
disk, which caused the bar to gain angular momentum. Such an event is 
rare; spirals generally remove angular momentum from the bar at most 
relative phases. It is possible that VK03 were unlucky to have such an 
event in their case, but they report similar behavior in their model B 
making a chance event unlikely. 
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One significant difference between our code and that used by VK03 
(Kravtsov, Klypin &; Khokhlov 1997) is that their resolution is adap- 
tive, which causes gravity to strengthen at short range when the grid 
is refined. The increase in the local density as the bar amplitude rises 
causes the code to refine the grid, strengthening gravity and thereby 
causing the bar to contract slightly and to spin-up. Wc have found that 
a reduction of softening length in our code at this epoch also leads to 
a metastable state. It is likely, therefore, that their anomalous result is 
an artifact of their adaptive code. 

6. The Metastable State is Fragile 

Whatever the origin of the bar speed-up in simulations, it remains 
possible that the metastable state could occur in real galaxies. If friction 
in a dense halo can be avoided for this reason, then the observed pattern 
speeds will provide no constraint on the halo density. 

However, further experiments in which we perturbed our model in the 
metastable state very slightly, revealed that the state is highly fragile. 
For example, a satellite of merely 1% of the galaxy mass flying by at 
30kpc is sufficient to jolt the system out of the metastable state. We 
therefore conclude that anomalously weak friction is unlikely to persist 
for long in nature. 

7. Conclusions 

Tremaine & Weinberg (1984) showed that angular momentum is trans- 
ferred from a rotating bar to the halo through resonant interactions. We 
find that friction is dominated by a single resonance at most times, and 
that corotation is most important for a bar with realistic pattern speed 
— i.e. when the bar extends almost to corotation. 

Friction arises because the phase space density is a decreasing function 
of angular momentum in normal circumstances, causing an excess of 
particles that gain angular momentum over those that lose. While this 
process would tend to flatten the density gradient if the pattern speed 
remained steady, the decreasing angular speed of the bar prevents this 
steady state from being reached. Instead we find that the density of 
particles in phase space develops a shoulder, with the resonance holding 
station on the high-angular momentum side of the shoulder as the feature 
moves to larger Lres- 

However, if the bar is spun up slightly for some reason after a period 
of normal friction, the rise in the pattern speed may move the reso- 
nance to the other side of the pre-constructed shoulder. The change in 
the local gradient of particle density at the dominant resonance causes 
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friction to become very weak for a while, allowing the bar to rotate al- 
most steadily. Mild friction persists because of contributions from other, 

sub-dominant resonances, and normal friction resumes once the pattern 
speed has declined sufficiently for the gradient at the main resonance to 
become favorable for friction once more. 

A state in which strong friction is suspended for this reason is "meta- 
stable", both because it relies on a local minimum in the phase space 
density, and because the state is fragile. A very mild jolt to the system 
is sufficient to cause normal friction to resume. 

The absence of friction in the simulation Ai reported by Valcnzucla 
&; Klypin (2003) is probably an artifact of their code. Their adaptive 
grid causes gravity to strengthen as the bar density builds up, mak- 
ing the pattern speed of the bar rise for a purely numerical reason. 
Thus their claimed counter-example to the argument of Debattista &: 
Sellwood (1998, 2000) is a numerical artifact of their method. Pace 
Holley-Bockelmann & Weinberg (2005), our constraint on halo density 
still stands: A strong bar in a dense halo will quickly become unaccept- 
ably slow through dynamical friction. 
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